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I. INTRODUCTION 



Microscopic descriptions of exotic nuclei are becoming more important in the light of 
recent experiments involving the scattering of radioactive ion beams (RIBs) off hydrogen. 
In inverse kinematics, such equate to scattering of protons from the ions, and with current 
methods of analysis [1] of such scattering, a most complete map of their matter densities 
can be made. Such well tested descriptions of exotic nuclei will be of relevance in analyses 
of data to be taken with proposed electron-ion collider being built at GSI [§], as well as from 
the SCRIT project at RIKEN [3J. Descriptions of the charge densities, especially of nuclei 
with neutron halos, will require detailed microscopic models to account for the structure of 
the core. 

Traditionally proton scattering has been one of, if not, the best means by which the matter 
densities of the nucleus may be studied. Microscopic models now exist that can predict 
results of both elastic and inelastic scattering reactions. When good, detailed specification 
of the nucleon structure of the nucleus is used, those predictions usually agree very well 
with observation, both in shape and magnitude. To facilitate such analyses of data, one first 
must specify the nucleon-nucleus (NA) interaction. To do so requires two main ingredients: 
a) an effective nucleon-nucleon (NN) interaction in-medium, allowing for the mean field 
as well as Pauli-blocking effects; and b) a credible model of structure for the nucleus that 
is nucleon-based. When the effective NN interaction is folded with the one-body density 
matrix elements (OBDME) of the target ground state, a microscopic NA interaction results. 
Such interactions have been used successfully in studies of the structures of stable nuclei 
(see |l| for a complete review of those studies), as well as of exotic nuclei j3, 111]. 

Herein, we consider the scattering of the exotic nuclei 8 He, 10 ' n C, and 18 > 20 > 22 O from hy- 
drogen, for which there are recent data. The data for the elastic scattering of 8 He from 
hydrogen at 15.7 A MeV [6[ was analyzed in terms of the JLM model using G matrix shell 
model wave functions of Navratil and Barrett Q, as well as a coupled-channels model in- 
volving coupling to the 8 He(p, cf) 7 He channel. Their analyses indicated a strong coupling 
to the (p, d) channel. This is problematic as 7 He is unbound. Any possible recoupling to 
reform 8 He in its ground state is highly unlikely. Indeed, their analysis based on this cou- 
pling required a scale strength of the imaginary part of the JLM potential of 0.2, which the 
authors attributed to compound nucleus effects. For 8 He, the energy (15.7 MeV) lies well in 
its continuum and compound nucleus contribution to the elastic scattering is not likely to 
be that large. So such a sizable reduction in the imaginary part of the JLM optical potential 
seems to be unrealistic. 

The JLM interaction was also used in the analyses of the scattering of the C [8| and O 
isotopes 0, [Io[ from hydrogen. For the analyses of the data from the C isotopes, and from 
i8,20q jlm potential was again used with the adjustment of parameters to fit the data. 
In the case of 22 [9j, two analyses were made: a) one using a folding model producing a 
real microscopic optical potential, but for which the imaginary part was obtained from a 
global phenomenological parameterization; and b) one using the JLM potential for which the 
neutron transition densities were adjusted to fit the data. As a consequence, those analyses 
may not necessarily be sensitive tests of the structures of the targets as any possible problems 
may be masked by a judicious use of normalisation factors. Thus we have reanalyzed those 
data seeking a better understanding of the structures of the exotic ions. To do so we have 
used the Melbourne g-folding model of scattering jl| together with a microscopic distorted 
wave approximation (DWA) to describe inelastic processes. The shell model has been used 
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to specify the putative structures of the ions. 



II. STRUCTURE AND REACTION THEORY 

In this section we give details of the structure of the exotic (RIB) nuclei and of the 
theories used to evaluate elastic and inelastic data from the scattering of those RIB nuclei 
from hydrogen. 



A. Many nucleon structures of He, ' C, and 



18,20,22 



o 



8 He lies at the neutron drip line. It is weakly bound with a two-neutron separation energy 
of 2.137 MeV It is also an example of a Borromean nucleus as 7 He is unbound, as is 
any two-body subsystem of 8 He. We have used a no-core shell model to define its ground 
state. That has been obtained in a complete (0 + 2 + 4)huj model space using the G-matrix 
interaction of Zheng et al. [12]. That ground state specification has been used previously in 
analyses of scattering data at 68 A MeV |4(. 

Spectra for 8 He are given in Fig. [TJ The (0 + 2 + 4)hu shell model results are compared 
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FIG. 1: The spectrum of 8 He. The result of the (0 + 2 + A)Hlo shell model calculation is compared 



with that of the VMC calculation 



not only with known data [IJ] but also with values determined by a variational Monte Carlo 



(VMC) calculation [14|. There is very little experimental information on the spectrum 
of 8 He. The first excited state is listed at 2.8 ± 0.4 MeV and has J W ;T = (2+);2 [ll|. 



Other states, at 1.3, 2.6, and 4.0 MeV [llj, have been suggested from a heavy ion transfer 
experiment, but as yet no other data support their existence. 

The results of the VMC calculation place the 2+ state in very close agreement with 
experiment. However, that calculation places an extra 1 + state in the spectrum. Not only 
has that state not been observed in experiment, but also it has not been found with other 
calculations. The shell model results clearly need improving to reduce the gap energy. 
Such suggest that an even larger space is required. However, binding energies reflect the 
asymptotic properties of wave functions and other tests are required to probe the credibility 
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FIG. 2: Energy spectrum of 10 C. The result of the (0 + 2)huj calculation is compared to the data 
of Ref. [lfi]. 



of current wave functions through the nuclear medium. A first test is the root mean square 
(r.m.s.) radii, which is most sensitive to characteristics of the (outer) surface of a nucleus. 
Using our shell model gave an r.m.s. value for 8 He of 2.63 fm, which is in good agreement 



with the value of 2.6 fm extracted from high energy data using a cluster model |15I . 

Wave functions of 10 C and 11 C were obtained in a complete (0 + 2)hu model space using 



the MK3W interaction (l6| . That model reproduced the 12 C spectrum to 20 MeV 17], with 
the exception of the superdeformed + state at 7.65 MeV. There is little known about the 
spectrum of 10 C as is shown in Fig. [2j The 2 + ; 1 excited state is at 3.354 MeV 18J while 



the SM prediction puts it at 4.272 MeV. The state at 6.58 MeV has been given a tentative 
assignment of 2 + for which the corresponding state in the SM spectrum lies at 5.613 MeV. 
There is no indication in the SM result of the two states below it. They have not been given 
assignments and may be intruder states. 

The spectrum of C is shown in Fig. [31 Therein, the results of (0 + 2)Huj and (1 + 3)/za> 
SM calculations, for negative and positive parity states respectively, are compared to the 
data [19|]. There is better agreement with the known spectrum in this case. 

The ground states for 18 > 20 > 22 O were obtained from a complete Ohu SM calculation using 
the USD interaction of Brown and Wildenthal 2(| . This model of structure suffices for use 
in describing the elastic scattering data from these nuclei. For the inelastic scattering to 
states in 20 O and 22 0, we also obtained the transition densities in this model, noting that 
there may be scaling of the cross sections due to core polarization corrections. 



B. The g- folding optical potential for elastic scattering 

A detailed study of the ^-folding optical potential has been published [l[ and so only 
salient features are presented. In coordinate space, that optical potential can be written as 

U(r,r';E) = 8(r-r')J2m [ ^«<? D (r, s; E)cp t (s) ds + ^ n^*(r)/(r, r'; E)^(r') (1) 



4 



E x (MeV) ? 

6 
5 
4 
3 
2 

1 





5/2 + 

7/2" '. 

3/2" 



5/2" 



1/2 



3/2"; 1/2 



expt. 



3/2 + 



MK3W 



FIG. 3: Energy spectrum of 11 C. The results of the shell model calculations are compared to the 
data of Ref. 191 ] . 



To evaluate these potentials requires specification of three quantities. They are the single 
nucleon bound state wave functions the orbit occupancies rti which, more properly 

are the nuclear state OBDME, and the NN ^-matrices g D ^ E (r, r'; E). Those g-matrices 
are appropriate combinations of NN interactions in the nuclear medium for diverse NN 
angular momentum channels [l|. For the latter, much success has been had using effective 
NN interactions, now commonly designated as the Melbourne force, which have the form 
= gfg (r, E\ kf(R)) where r = |ro — r^J and R = | (r + r^). Therein the Fermi momenta 
relate to the local density in the nucleus at distance R from the center when are the 
coordinates of the colliding projectile and bound nucleons. {ST} are the iViV spin and 
isospin quantum numbers. 

For use in the DWBA98 program, these effective NN ^-matrices are, more specifically, 
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(2) 



where 0j are 
{1, (<t-<t),(t-t), , 
body spin-orbit force 



[a 



the characteristic operators for central forces (i = 1), 
or ■ t)}, for the tensor force {% = 2), {Si 2 }, and for the two- 
= 3), {L • S}. The Sj (E;kf) are complex, 



energy- 



and 



density-dependent strengths. The properties of the ^-matrices are such that the ranges of 
the Yukawa form factors, as assumed in the above, can be taken as independent of energy 
and density [l|. 

The strengths (and ranges) in these effective iViV interactions were found by mapping 
the double Bessel transforms of them to the iViV ^-matrices in infinite nuclear matter that 
are solutions of the partial wave Bethe-Brueckner-Goldstone (BBG) equations [l|. With the 
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BBG g-matrices denoted by #bbg;LL'(9'' 1' E), the mapping is 

poo 

9Z$G;lM, 95 E) = Vi) / ^ X 3v (q'r) <$ OT (r, fc,) j L for) 

(>) 

= J2 ft) s f M / ^ 2+ V tfO ^ (9r) ^ 

V ./o 



(3) 



where a : {LL'JST} and A = save for the tensor force for which it is 2. In application, a 
singular-valued decomposition has been used to effect this mapping and it was found that 
four Yukawa functions for the central force and four with other ranges for the spin-orbit and 
tensor forces sufficed to give close mapping to both on- and, a near range of, off-shell values 
of the BBG (/-matrices in 32 iViV S, T channels and for energies to 300 MeV. 

The other requirements come from the assumed models of structure. The many-nucleon 
aspects, the OBDME, are the ground state reduced matrix elements of particle-hole opera- 
tors 



S j <; 



g.s. 



a\- x a,j 



g.s. 



(4) 



and which are defined more generally in regard to transitions in the following subsection. If 
there are no non-Hartree-Fock contributions, these OBDME reduce to the shell occupancies 
in the ground state. The single nucleon (bound state) wave functions usually are chosen 
either to be harmonic oscillators with oscillator length as used in the shell model calculations 
or from a Wood-Saxon potential. The parameters of that Wood-Saxon potential are chosen 
to give wave function solutions that meet some criteria such as the r.m.s. radius or an 
electron scattering form factor. 



C. The distorted wave approximation for inelastic scattering 

In the DWA, amplitudes for inelastic scattering of nucleons through a scattering angle 8 
and between states Jj, Mj and Jf, Mf in a nucleus, are given by 

1 J f J i \P) = 

X { j\k o 0) (* JfMf (l ■■■A)\AJ2 9 S Jir QU E- kf) P s P T A 01 {\xi + \m) 1*^(1 -"A))}, 

(5) 



ST 



where u, v' are the spin quantum numbers of the nucleon in the continuum, x^ are the 
distorted waves, and g^ff (roi, E; kf) PsPt is the spin-isospin Melbourne force. The operator 
Am effects antisymmetrization of the two-nucleon product states. 
Then, by using cofactor expansions of the nuclear states, i.e. 

\®Jm) = ~Ff5Z IVW 1 )) a jm \VjAf) , (6) 
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the matrix elements become 



rpMfMiW 1 

j f Ji 



The density matrix elements in the amplitudes reduce as 



(0) = E {JfMf a) 2m2 a nmi 
hhST 

x ( x { j\K0)\ & j2m2 (l)\g*I (r„i, E- k f ) PsPtAoi {\xl +) (M) |^- imi (l)>} 



(7) 
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(8) 



where Sj ± j 2 j are the transition OBDME. The DWA amplitudes are then (with £ 
{ji, 32,m u m 2 ,I(N),S,T}) 



2J f + 1 



Siiia/ 0i h ~m 2 \IN) (Ji Mi I N\ J f M f ) 



g v ~"f 

(-) 



Xi7 (M) (^(l)|flS'(ro 1 ,^;fc/)P5^^{|xS f) (ki0)>|^ imi (l)>} . (9) 

In our calculations of these DWA amplitudes, a) the g-folding model has been used to 
determine the distorted waves in both the incident and emergent channels, b) the same 
effective iViV interaction has been used as the transition operator, c) shell models have been 
used to find the transition OBDME, and d) the single nucleon wave functions chosen to 
form the optical potentials have also been used for the single nucleon states. 



III. RESULTS 



In previous papers [H|, [2l| , the cross sections from the elastic and inelastic scattering (to 
the 2f state) of 40. 9A and 24.5A MeV 6 He ions from hydrogen were analyzed using the 
methods described above. Those results provided clear evidence that 6 He had an extended 
neutron matter distribution consistent with a neutron halo. Similar analyses of elastic 
scattering data from the scattering of 72A MeV 8 He ions from hydrogen suggested that it 
had a neutron skin but not the extended distribution of a halo |4J. Recently [6| data from 
the scattering of 15. 7A MeV 8 He ions from hydrogen has been reported, and we consider 
that now. The energy is the smallest at which the ^-folding method has been applied. In 
most nuclei in the excitation energy range ~ 10 to ~ 30 MeV there are giant resonances; the 
coupling to which is known (22j to influence nucleon scattering. When that is so, corrections 
to the (^-folding premise need be made. However, an important exception is the set of 
Helium isotopes for which there are no giant resonances. With 8 He also, the nucleon break- 
up threshold lies at just a few MeV so that 15 MeV is well in the continuum of the nucleus. 
For these reasons, it should be appropriate to use the (7-folding approach to analyze the 
15. 7 A MeV data. 
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A. The elastic scattering of 15.7^4 MeV He from hydrogen 

Using the structure of the ground state described previously and with the 15 MeV Mel- 
bourne force, two predictions of the 15. 7 A MeV 8 He-p scattering have been made. The 
results are compared with the data in Fig. |H The result shown by the dashed curve was 
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FIG. 4: (Color online) Differential cross sections for the elastic scattering of 15. 7 A MeV 8 He ions 
from hydrogen. The data 0] are compared with g-folding model results. The curves are identified 
in the text. 

found by using oscillator wave functions for the single nucleon states in the folding process, 
while that shown by the solid curve was found using Woods-Saxon functions generated with 
the potential used before [I, Clearly the predictions are in very good agreement with 
this data with preference (as was found before) for the Woods-Saxon single-nucleon wave 
function set. The JLM model results [o| do better in reproducing this data around 80° 
but do not as well at forward scattering angles. However, the JLM results required some 
a posteriori adjustment of the real and imaginary strengths of the optical potential from 
those considered standard. That, and the strength of the measured 8 He(p, d) cross section, 
led them to believe that a second order (p, d) (d, p) process was important in the description 
of elastic scattering. But, the measurement of the (p, d) reaction was inclusive with no iden- 
tification of the final state being 7 He. As 7 He is unbound against neutron emission there 
are many multi-particle exit channels from the (p, d) reaction and the likelihood of the con- 
glomerate reforming to a 8 He and proton is very small. In fact the plethora of propagating 
channels (given the energy and momentum sharing possible) is one reason why an optical 
potential approach is appropriate. 

B. The scattering of 1011 C ions from hydrogen 

Another GANIL experiment [8J] determined differential cross sections for the elastic and 
inelastic scattering of 45. 3A MeV 10 C ions and of 40. 6A MeV n C ions from hydrogen. 
That data are compared with predictions in Figs. [5] and [61 Previously, data taken at 40. 9A 
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MeV [5|, Hll of 6 He ions elastically and inelastically scattered from hydrogen were predicted 
well with the g-folding method. The current results, with which the new data are compared 
in these two figures, were found by using the structure as described earlier. Also we have 



used the same Melbourne force that was used in the analyses of the 6 He data 2l[ and 
either oscillator or Woods-Saxon wave functions for the bound nucleons. Note that in these 
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FIG. 5: (Color online) Differential cross sections for the elastic (top) and inelastic scattering to the 
state (bottom) of 45.3A MeV 10 C ions from hydrogen. 




FIG. 6: (Color online) Differential cross sections for the elastic (top) and inelastic scattering to the 
§~ state (bottom) of 40.6A MeV n C ions from hydrogen. 

diagrams the inelastic scattering data and results are plotted on a linear scale to emphasize 
any shortfalls. 
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In Fig. [5j cross sections for the scattering of 10 C are shown. Clearly the elastic scattering 
prediction agrees with the data almost as well as the phenomenological JLM analyses of 
Jouanne et al. |8|. With the inelastic scattering (to the 2 + ; 3.36 MeV state), the results found 
using the OBDME from the shell model calculations described before and with oscillator 
(dashed) and Woods-Saxon (dot-dashed) functions for the single-nucleon states required 
in the DWA calculations lie below the data values. A scaling of 1.5 on the Woods-Saxon 
result gives the solid curve which is in very good agreement with the data. Scaling was 
also required with the phenomenological model analyses In our case, however, given the 
propriety (admittedly mostly justified through numerous uses) of 

1. the model, and in particular with it, specific treatment of the exchange amplitudes; 

2. the transition operator (Melbourne force) in which account is taken of medium effects 
on the NN interaction; and 

3. the single-nucleon wave functions which link well to ground state properties; 

the shortfall we attribute to a limitation of the many-nucleon structure for the ion itself. It 
has been long known [1] that core polarization corrections, be they made phenomenologically 
or by appropriate increase in the space in which the structure is assessed, influence inelastic 
scattering evaluations strongly; and most obviously when collective enhancing excitations 
are involved. That the amount in the case of 10 C is but a 50% enlargement attests to the 
reasonable first guess that is given by the structure calculations made using the MK3W 
potentials. The enlargement equates to a polarization charge, Ae ~ 0.1 le. 

In Fig. El the 40. 6A MeV n C results are shown. The elastic scattering data are again well 
predicted by the ^-folding method. Likewise the cross section for excitation of the | ; 4.32 
MeV state as predicted using Woods-Saxon wave functions (dot-dashed curve) reproduces 
the observed data very well when a small upward scale (of 1.4) is made. 

Thus the MK3W structure adopted for these states in 10 ' 11 C is quite good. Use of that 
structure in the scattering calculations gave cross-section shapes very like that observed in 
the data and requiring but a small amount of core polarization to reproduce the measured 
data from the inelastic (E2 dominant) transitions. We surmise that the first complete basis 
space enlargement possible to the shell model scheme we have used to describe these nuclei 
will give much, if not all, of such enhancement. 



C. The 18-20,22q isot0 p es 

Recently the cross sections from 43A MeV 18 > 20 O ions and from 46. 6A MeV 22 ions 
scattering from hydrogen targets have been measured 0, [l(| • Cross sections for both the 
elastic and the inelastic excitation of 2f states of those ions were obtained. Those data are 
shown in the top, middle, and bottom panels of Fig. [7J Those data are compared with our 
predictions (^-folding for the elastic and DWA for the inelastic transitions) built with the 
Ohio shell model wave functions and the Melbourne force. The elastic scattering cross-section 
predictions depicted by the solid lines agree well with the data (filled circles) save that the 
18,20(3 results have an over-pronounced minimum near 40°. The inelastic scattering cross 
sections given by our approach are those depicted by the long dashed curves and they all 
are smaller than any of the data sets (opaque circles). The solid curves that match the 
inelastic data well are those calculated values increase by factors of 2 ( 18 0), 5 ( 20 O), and 
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FIG. 7: (Color online) Differential cross sections for the elastic and inelastic scattering (to the 
states) with the isotopes of oxygen, top ( 18 0), middle ( 20 O), and bottom ( 22 0). The energies are 
as indicated. 



1.6 ( 22 0). These ratios are quite similar to those of the squares of the deformations required 
in a previous analysis of this data 23|. That analysis also found very good fits to all of 
the elastic scattering data, but to do so involved using many parameters including separate 
renormalization of the real and imaginary components of the optical potentials, and for each 
nucleus independently. 

With our results for the inelastic scattering cross sections one can envisage different 
degrees of core polarization being required with the base valence models having 2, 4, and 6 
(full sub-shell) (is -orbit occupancies as the dominant term in the description of the ground 

states of 18 ' 20 ' 22 O respectively. That 20 O requires largest scaling is reflective of that nucleus 
having considerably more configuration mixing than the other isotopes in a good description 
of its states. 



IV. CONCLUSIONS 



The microscopic ^-folding approach to analysis of data from the elastic scattering of ra- 
dioactive ion beams from hydrogen targets when the beam energy equates to an excitation 
of those ions sufficiently in their continua, is most appropriate. When essentials are satis- 
fied, that approach can, and does, give predictions of elastic scattering that one may have 
confidence in being a good representation of observation. It is essential, though, to use a 
realistic energy- and medium-dependent effective NN interaction in the folding. Also it is 
essential that as good a prescription as possible be made for the distributions of nucleons in 
the ground state of the ion. Finally one must be faithful to the Pauli principle, not only in 
the specification of the effective NN interaction, but also in the calculation of the scattering 
allowing for the knock-out process. Doing so introduces a specific non-local term in the 
optical potential. Our experience is that such non-locality should be treated as exactly as 
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possible and not localized to simplify the problem to solution of a set of local Schrodinger 
equations. 

We have also shown that the DWA suffices to analyze the inelastic scattering cross sec- 
tions. There is a proviso that the energy is sufficient that coupled-channel effects due to 
discrete and/or giant resonance states of the ion are not important. When such is not the 
case, other methods of analysis, such as with a multichannel algebraic scattering theory [24] , 
are more relevant. Even then it has been noted that due care of the Pauli principle is es- 
sential 25|]. But for data at the energies considered, the DWA worked very well. With all 
details preset, the DWA results are also predictions which, in the cases considered, well re- 
produced the shapes of the data. All of those calculated results had to be enhanced to meet 
the magnitudes of that data but those scalings were consistent with the core polarization 
one must expect for shell model wave functions defined in the small bases we have used. 
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